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ON THE RATE OF CONVERGENCE OF SIMPLE AND 
JUMP-ADAPTED WEAK EULER SCHEMES FOR LEVY 

DRIVEN SDES 

R. MIKULEVICIUS 



Abstract. The paper studies the rate of convergence of a weak Eu- 
ler approximation for solutions to possibly completely degenerate SDEs 
driven by Levy processes, with Holder-continuous coefficients, ft inves- 
tigates the dependence of the rate on the regularity of coefficients and 
driving processes and its robustness to the approximation of the incre- 
ments of the driving process. A convergence rate is derived for some 
approximate jump-adpted Euler scheme as well. 



1. Introduction 

The paper studies the weak Euler approximation for solutions to possibly 
completely degenerate SDEs driven by Levy processes. As in [12] . the main 
goal is to investigate the dependence of the convergence rate on the regularity 
of coefficients and driving processes. In addition, we consider the robustness 
of the results to the approximation of the law of the increments of the driving 
noise in the whole scale of time discretization errors. 

Let (f2, J 7 , P) be a complete probability space with a filtration F = {^i}te[o,Tl 
of cr-algebras satisfying the usual conditions and a G (0, 2] be fixed. Consider 
the following model in R d : 

(1.1) X t =X + [ a(X s )ds+ [ b(X s )dW s + [ G(X s -)dZ s ,t G [0,T], 
Jo Jo Jo 

where a(x) = (a l (a;))i<i<d> H x ) = ( x ))i<i<d,i<j<n, G(x) = (G % i (x))x<i<d,i<j< m , 
x G R d are measurable and bounded, with a = if a G (0, 1) and b = if 
a G (0,2). The process W s is a standard Wiener in R". The last term is 
driven by Z = {Zt\t&\o,Th an w-dimensional Levy process whose character- 
istic function is exp{tr/(^)} with 
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where x a ( v ) = X{|«|<i} 1 {ae(l,2]}- Hence, 



Zt = J J {I ~ x a ( v )) v P(ds,dv) + J J x a ( v ) v <l(ds,dv), 

where p(dt, dv) is a Poisson point measure on [0, oo) x R™ (R™ = R m \{0}) 
with E\p(dt,dv)] = n(dv)dt, and q(dt,dv) = p(dt,dv) — n(dv)dt is the cen- 
tered Poisson measure. It is assumed that Z t is a Levy process of order 
a : 



(\v\ a A l)n{dv) < oo. 

Let the time discretization {rj,z = 0, . . . , n^} of the interval [0, T] with 
maximum step size 5 > be a partition of [0, T] such that = tq < t\ < 
• • • < T nT = T and maxj(rj — Tj_i) < 5. The Euler approximation of X is 
an F-adapted stochastic process Y = {^t}tg[o,T] defined by the stochastic 
equation 

(1.2) Y t = X + I ' a(Y Tis )ds+ I b(Y T JdW s + f G(Y T JdZ s ,t e [0,T], 
Jo Jo Jo 

where Tj s = Tj if s € [t^, Tj+i), i = 0, ...,tit — 1. Contrary to those in 
(jl.ip . the coefficients in (|1.2p are piecewise constants in each time interval 

Of [Ti,T i+1 ). 

The weak Euler approximation Y is said to converge with order k > if 
for each bounded smooth function g with bounded derivatives, there exists 
a constant C, depending only on g, such that 

\Bg(Y T )-Eg(X T )\ <C5 K , 

where 5 > is the maximum step size of the time discretization. 

The weak Euler approximation of stochastic differential equations with 
smooth coefficients and G = has been consistently studied. For diffu- 
sion processes, Milstein was one of the first to investigate the order of weak 
convergence and derived k = 1 [131 E] . Talay considered a class of the sec- 
ond order approximations for diffusion processes [El H9] . For Ito processes 
with jump components (a finite number of jumps in a finite interval), it was 
shown in !)j the first-order convergence in the case in which the coefficient 
functions possess fourth-order continuous derivatives . Platen and Kloeden 
h Platen studied not only Euler but also higher order approximations [5l PT5] 
and references therein. 

Protter &: Talay (|17j) analyzed the weak Euler approximation for (|l.ip 
with a = 2. They proved that the order of convergence is k = 1, provided 
that G, b, a and g have four bounded derivatives and the Levy measure of 
Z has finite moments of the order fi = 8. In this paper we show that k = 1 
can be achieved when fi = 4 and there still is some order of convergence for 
fi € (2,4]. Moreover, we assume /3-Lipshitz continuity of the coefficients and 
g and derive that for a < f3 < fi < 2a the order of convergence k = ^ — 1. In 
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particular, when /3 = fj, = 2a with a € (0,2) (the diffusion part is absent), 
the convergence order is still k = 1. 

As in [10] and [12], this paper employs the idea of Talay (see [JB]) and 
uses the solution to the backward Kolmogorov equation associated with Xt, 
Ito's formula, and one-step estimates. Since one step estimates were derived 
in |12j . the main difficulty is to solve the degenerate backward Kolmogorov 
equation in Lipschitz classes (see Theorem |4] below). We obtain the solution 
of the degenerate equation as a limit of solutions to regularized (nondegen- 
erate) equations. Although the solution to (11.11) is strong and probabilistic 
arguments are applied for the uniform Lipshitz estimates of the approximat- 
ing sequence, contrary to [T7], we do not use derivatives of the stochastic 
flows. 

If (jl.ip has a nondegenerate main part, some assumptions imposed can 
be relaxed (see [12], |10| . Kubilius & Platen and Platen & Bruti-Liberati 
[HQS]). More complex and higher order schemes were studied and discussed, 
for example, by Cont and Tankov, Jourdain and Kohatsu-Higa (see [JJ, [3] 
and references therein). 

Motivated by the difficulty to approximate the increments of the driving 
processes, Jacod, Kurtz, Meleard and Protter in |3J, studied the approx- 
imated Euler scheme where the increments of Z are substituted by i.i.d. 
random variables that are easier to simulate. There are two sources of er- 
rors in this case. One comes from time discretization and the other one from 
substitution. We extend some of the results in [3] to the whole rate scale 
and show that the errors add up. In particular, the driving process Z can 
be replaced with a Levy process Z having finite number of jumps in [0, T] 
by possibly cutting small jumps of Z and sometimes replacing them with 
a Wiener process or drift. In addition, we consider a simple jump-adapted 
Euler scheme and show that presence of Z-jump moments in the partition 
{Tj} influences the convergence rate. The approximation itself is simpler and 
assumptions imposed are different than those introduced by Kohatsu-Higa 
and Tankov in [6] (see references therein as well) for a more sophisticated 
(higher order) jump-adapted scheme. 

The paper is organized as follows. In Section 2, some notation is intro- 
duced, the main results stated and the proof of the main theorem is outlined. 
In Section 3, we present the essential technical results about backward de- 
generate Kolmogorov equation, followed by the proof of the main theorem 
in Section 4. The robustness of the approximation and jump-adapted Euler 
scheme is considered as well. In the last section we discuss the optimality 
of the imposed assumptions. 



2. Notation and Main Result 

Denote H = [0,T] x R d , N = {0,1,2,...}, R$ = R d \{0}. For x,y e H d , 
write (x,y) = J2i=i x iUi- F° r (*> x ) £ H, multiindex 7 £ N d with D 1 = 
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d 



Au(t,x) = y~]u XiXi (t, 



dx ji dx id , and i,j = l,...,d, denote 
d 

d t u{t, x) = —u(t, x), D k u(t, x) = (D^u(t, x)) |7|=fc , k G N, 

diu(t,x) = u Xi (t,x) = -Q—u(t,x), d$ju(t,x) = u XiXj (t,x) = u(t,x) 
d x u(t,x) = Vu(t, x) = V x u(t, x) = [d\u(t, x), . . . , d^u{t, x)) , 

d 

i=i 

For a smooth function v on R d and /c G N, denote 

d 

v^( X -,e,---,e)= e «* il ..^ fc (x)4...et a; ^ <GR<i » i = i .-' fc - 

ii,...,i fc =l 

In particular, £) = (Vv(x),^),x, £ € R d . 

For /3 = [/?]- + {/?}+ > 0, where [0]" G N and {/?}+ G (0, 1], let C^(iT) 
denote the Lipschitz space of measurable functions u on H such that the 
norm 

|u|^= ^ |£>M*,*)lo + sup 

i7i<[/3i- w=gr. 

is finite, where |u|o = sup/ ta ,) 6i? We denote C^(R rf ) the correspond- 

ing function space on R d . 

C = C(-, . . . , •) denotes constants depending only on quantities appearing 
in parentheses. In a given context the same letter is (generally) used to 
denote different constants depending on the same set of arguments. 

The main result of this paper is the following statement. 

Theorem 1. Let a < j3 < \i < 2a, and assume a\b ij G C^(R d ), G ij G 
C^ vl (R d ), and 



\v\ a dTT + / \v\L l ir(dv) < 00, 
M<i J\v\>i 

where ir is the Levy measure of the driving process Z. Then there is a 
constant C such that for all g G C^(R d ) 

\Eg(Y T ) -Bg(X T )\ < C\g\^-\ 

Applying Theorem [1] to the case a = 2 we have an obvious consequence 
in the jump-diffusion case. 

Corollary 1. Consider the jump- diffusion case (a = 2) 

X t = X + [ a(X s )ds+ I b(X s )dW s + [ G(X s _)dZ s ,t G [0,T]. 
Jo Jo Jo 
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Let 2 < < [i < 4. Assume a, b ij , G ij G C p (R d ) and 



\v\ ir(dv) + / \v\ lx Tr(dv) < oo. 

\v\<l J\v\>l 

Then there is a constant C such that for all g G (7^(R d ) 

\Eg(Y T ) -Bg(X T )\ < C\g\^-\ 

An immediate extension of Theorem [I] (for the test function g G (^(R ) 
with v G (0,/?]) is the following statement. 

Corollary 2. Let a < f3 < ^ < 2a, and assume a\b ij G CP(R d ), G ij G 
C^ vl (R d ), and 



/ \v\^Tt(dv) < oo, 

J|u|>l 



\v\ a dTT + 
'\v\<l J\v\l 

where ir is the Levy measure of the driving process Z. Let v G (0, f}\. Then 
there is a constant C such that for all g G C u (R d ) 

\Eg(Y T ) -Bg(X T )\ < C\g\ v 8 v( ^ ] ■ 

Remark 1. Ln particular, if a G [1,2], fi = f3 = 2a and g is Lipshitz 
(v = 1), then the convergence rate k = ^j. 

2.1. Approximate simple Euler scheme. Following [3j, for a G (0, 1), 5 > 

0, we choose a time discretization {-Tj} and replace the increments of the 
driving process Z Ti+1 — Z Ti in ll.2l bv T Ti -conditionally independent random 
variables Q, i = 0, . . . , tit — 1- We assume that there is a function cp(a) such 
that limo-^o 4>( a ) = an d for i = 0, . . . , nr — 1, 
(2.1) 

\B[h(Z Ti+1 - Z Ti ) - h(( i+1 )\T Ti }\ < C\h\ p <K*)(T i+ i - n), h G &(R d ). 

with some constant C, independent of a, 5 and h. Let £ t = if < t < 
T i-> = d ^ U < t < ti+iii = 1) ■ ■ ■ j ray— 1. We still assume that maxj(rj + i — 
T~i) < 5 and approximate X t by 

(2.2) f t = X + /' a(f Tis )ds + T 6(f Tis )aW s + / G(Y Tis )d£ a ,t G [0, T]. 
Jo J JO 

In this case 1^ depends on 5 and <r. 

In the following example we approximate the increments of Zt by the 
increments of a Levy process with finite number of jumps in [0, T]. This 
approximation is constructed by cutting small jumps of Z%. We replace 
the small jumps part by appropriately chosen drift if a < /3 G (1, 2], a G 
(0,1]. If a < f3 G (2, 3], a G (1,2], the small jumps part is replaced by 
a Wiener process. Given a G (0,1), we denote B a the square root of the 



positive definite mxm -matrix ( J, ViVjdn I . Let Wt be a standard 

V I I— CT / l<i,j<m 

independent Wiener process in R m . 
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Example 1. For a £ (0,1) we approximate 

z t = J j "(1 - X a (v))vp(ds,dv) + J J x a (v)vq(ds,dv),t G [0,T], 

by 

z t = Z>1 + R?t 

with 

Zt = [ [ (l-x a (v))vp(ds,dv)+ [ [ x a (v)vq(ds,oIv) 

JO J\v\>a JO J\v\>a 

and 

( tf M<a v7r(dv) ifa<P& (1,2], a G (0,1], 



Rt= < 



l\v\<_ 

B°W t i/a</3G(2,4],aG(l,2], 



otherwise. 
In this case (see Lemma\E below) \2. 1\) holds with 

<P{a)= f \vf A3 dn 

J \v\ <<J 

and 

(2.3) Ci+i = Z Tl+1 -Z Ti ,i = 0,...,n T -l. 

We show that time discretization and substitution errors add up. 

Theorem 2. Let a < f3 < fj, < 2a, and let a\W G &(R d ), G ij G 
C^ vl (R d ), and 



/ \v\ a dir + / \vY L it(dv) < oo, 

J\v\<l J\v\>l 



'I 

where tt is the Levy measure of the driving process Z. Assume that there is 
a function <ft(o~) such that lim^^o 0( cr ) = an d for i = 0, . . . , ut — 1, 

(2.4) |E/i(Z Ti+1 -Z Ti )-E/»(C i+1 )| <C\h\ /3 J>(o-)(T i+1 -T i ),heCP(R d ), 

for some constant C . 

Then there is a constant C (independent of a, 5) such that for all g G 

&{n d ) 

|E<7(F T ) - E 5 (X r )| < Cl^i^^i- 1 + </>(?)]. 

The same way as Corollary [2] (see the proof below) we have the following 
statement. 

Corollary 3. Let assumptions of Theorem^ hold and v G (0,/?]. Then 
there is a constant C such that for all g G C u (R d ) 

\Eg(Y T ) - Bg(X T )\ < C\g\ v [S^«-^ + 0(a)?]. 



CONVERGENCE RATE OF WEAK EULER SCHEMES 7 

Remark 2. (i) Assume the assumptions of Theorem\^hold. Since linio-^o i'i ') 
0, for each 5 > there is a = a (5) such that (f>(a(5)) < and therefore 

\Bg(Y T )-Bg(X T )\ <C\g\^-\ 

In particular, if 4>(o~) < Ca^ with fj, > (it is the case in Example U\ for 
a small jumps a' -stable-like driving process Z with a' < a), then we can 

choose cr M = <5q _1 or a = S^ -1 ^ . 

(ii) In order to study precisely the case of unbounded test functions (like 
one in one would need to solve first the backward Kolmogorov equation 
in Holder spaces with weights that are defined by the powers of w{x) = 

(l + |x| 2 ) 1/2 ,xGR d . 

Applying Theorem [2] to the model of Example [I] we have 

Proposition 1. Let a < /3 < // < 2a, and let a l ,b ij G &(R d ), G ij G 
C^ vl (K d ), and 



v\ a dTT + / \v\ fl TT(dv) < OO, 

v\<l J\v\>l 

where ir is the Levy measure of the driving process Z. For the approximate 
Euler scheme in Example^ there is a constant C (independent of a, 5) such 
that for all g G C^(R d ) 

\Bg(Y T ) - Bg(X T )\ < C\g\^~ l + ! \vf A3 d7v]. 

J \v\<a 

2.2. Approximate jump-adapted Euler scheme. As in Example HJ for 
a G (0, 1) we approximate the increments of the driving process 

Z t = j J (1 - x a (v))vp(ds,dv) + J J x a (v)vq(ds,dv),t G [0,T], 
by the increments of 

Zt = %t + -Rf> 

with 

Z? = f [ (1 - x a (v))vp(ds,dv) + [ [ x a (v)vq(ds,dv) 

JO J\v\>a JO J\v\>a 

and 

{ t !\v\<a v ^ dv ) ifa</3G(l,2],«G(0,l], 

Rt = { B a W t if a < /3 G (2,4], a G (1,2] 

otherwise, 
where B a is the square root of the positive definite m x m -matrix (^j^ <cr ViVj 
and Wt is a standard independent Wiener process in R r ' 



dir 
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Given a £ (0,1), 5 > 0, consider the following Z°"-jump adapted time 
discretization (see [9]): to = 0, 

(2.5) T i+ i = inf (t > n : AZ? + 0) A (t 4 + 5) A T. 

In this case the time discretization {Tj,i = 0, . . . ,ny} of the interval [0,T] 
is random, Tj are stopping times. We approximate X t by 

(2.6) Y t = X + f a {Y Tis )ds+ f* b(Y T JdW s + f G(Y T JdZ s ,t e [0,T]. 

Jo jo Jo 

The following error estimate holds. 

Theorem 3. Let a < p < fi < 2a, and let a\b ij € C p (R d ), G ij € 
C^ 1 (R d ), and 



\v\ a dir + / \v\^-K{dv) < oo, 
v\<l J\v\>i 

where ir is the Levy measure of the driving process Z. 

Then there is a constant C (independent of a, 5) such that for all g € 

\Eg(Y T ) - Eg(X T )\ < C\g\ p [{(S A A^ 1 )^}^ 1 + / \vf A3 d7r], 

J\v\<a 

where A CT = it ({\v\ > a}) and 

X a = 1 + l ae (i j2 )| / vdn\. 

Ja<\v\<l 

In particular, the following statement holds. 

Corollary 4. Suppose the assumptions of Theorem{3\hold. 

(i) Lf 5 = T (only jump moments are chosen for the time discretization), 
then 

\Bg(Y T ) - Eg(X T )\ < C\g\J ^ ) I" 1 + / \vf A3 dn]. 

\-V J J\v\<cr 

(ii) //sup CTe(01) | X T< j ti j< 1 vdir\ < oo for a € (1,2), then 

\Eg(Y T ) - Eg(X T )\ < C\g\p[(6 A A" 1 )" -1 + / \vf A3 dir], 

J\v\<(7 



where X a = ir ({\v\ > a}) . 



2.3. Outline of Proof of Theorem Q3 To prove Theorem[TJ as in [10] and 
|12j , the solution to the backward Kolmogorov equation associated with Xt is 
used. First we introduce the operator of the Kolmogorov equation associated 
with Xt- 
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For u G &P(H),P > a, denote 

1 d 

L z u{t,x) = (a(z),V x u(t,x)) + -J2(b\z),^(z))df j u(x) 

[u(t,x + G(z)v) - u(t,x) - x a ( v )(yxu(t,x),G(z)v))]ir(dv), 

Lu(t,x) = L x u(t,x) = L z u(t,x)\ z=x , 
where b % (z) = (b %j (z))i<j< m , i = l,...,d. 

Remark 3. Under assumptions of Theorem^ there exists a unique strong 
solution to equation (11.11) and the stochastic process 



u{X t )- f Lu(X s )ds,VueCP(R d ) 
Jo 



with f3 > a is a martingale. The operator L is the generator of Xt defined 
in (fTTT) . 

If v(t,x), (t,x) G H satisfies the backward Kolmogorov equation 
(d t + L)v(t,x) = 0, 0<t<T, 
v(T,x) = g(x), 

then by Ito's formula 

E[g(Y T )]-E[g(X T )} = E[v(T,Y T )-v(0,Y )]=B[ F (d t +Ly )v(s,Y s )ds\. 

Jo 

The regularity of v determines the one-step estimate and the rate of conver- 
gence of the approximation. 

3. Backward Kolmogorov Equation 

In Lipshitz spaces C^(H), consider the backward Kolmogorov equation 
associated with Xt. 

(3.1) (d t + L)u(t,x) = f(t,x), 

u{T,x) = g{x). 

Definition 1. Let f,g be measurable and bounded functions. We say that 
u G (JP{H) with ft > a is a solution to (|3.ip if 

(3.2) u(t, x) = g(x) + J [Lu(s,x) - f(s,x)]ds,V(t,x) G H. 

First we show that L:C"(H) — > C^~ a {H) is continuous. 
Lemma 1. Let a < j3 < \x < 2a, 

\v\ a dir + / {v^dir < oo 
v\<l J\v\>l 
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and a\W G &(R d ), G ij G C^ x {R d ). Then for any v G C p (R d ) we have 
Lv G CP- a (R d ) and there is a constant independent of v such that 

\Lv\/3- a < C\v\p. 

Proof. Let 

Bv(x) = / [v(x + G (x) v ) — v (x) — Xa( v )(^ v ( x )i G(x)v )] dir. 



Then 

Lv = Bv + (a{x),Vv(x)) + -(b\x), b> (x))d 2 Mx) . 

2 J 

By Proposition 13 in [[12], Bv G C^- Q (R d ) if /3-a £ N and \Bv\p_ a < C\v\p. 
In this case, obviously, Lv G C^~ a (R d ) as well. 
If a > = 1 + a, then 

Bv(x) = f [ [Vv(x + sG(x)v) - Vv(x)]G(x)v]dsdir 

J\v\<lJo 



Since 



+ / [v(x + G(x)v) — v(x)] dir. 

'\v\>l 



V(Bv(x)) = [ [ [d 2 v(x + sG(x)v) - d 2 v(x)]G(x)v}dsdTr 

J\v\<l JO 

+ [ [ d 2 v(x + sG(x)v)VG(x)v G(x)vdsdn 
J\v\<l JO 



+ / [Vv(x + G(x)v) -Vv(x)]dTT 

J\v\>l 

+ / Vv(x + G(x)v)S7G(x)vdir, 
J\v\>l 



>\v\>l 

it follows that sup x \V(Bv(x))\ < C\v\p. Therefore \Lv\p_ a < C\v\p as well. 
If a = 1 and ft = 2, then 

|V5u(x)| = J[Vv(x + G(x)v) -Vv(x)]dir + J Vv(x + G(x)v)G(x)vdir, 
sup\VBv(x)\ < C\v\ p 



X 



and I-LdI^q, < Cj-ul^. The case j3 = 4, a = 2 is considered in a similar 
way. □ 

The main result of this section is the following statement. 

Theorem 4. Let a < /3 < \x < 2a, and 



\v\ a ir(dv) + / |u| M 7r(df) < oo. 
u|<l J\v\>l 
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Assume a\W G C^(R d ),G^ G C^ vl (R d ). Then for each f G &{K d ),g G 
C^(R d ) ; there exists a unique solution u G C^{H) to (|3.ip and a constant 
C independent of f,g such that \u\p < C(\f\p + \g\p). 

To prove Theorem for e G (0, 1) we consider a non-degenerate equation 

(3.3) (d t + L e )u(t,x) = f(t,x), 

u(T,x) = g £ (x), 

where L £ u = —e a (— A) Q//2 u + Lu and 

9e{x) = / g(y)w £ (x - y)dy = / g(x - y)w e (y)dy,x G R d 



with = £~ d ™ (|) , x G R d , w G C^(R d ), f wdx = 1. 

An obvious consequence of Corollary 9 in [12] is the following statement. 



Lemma 2. fsee Corollary 9 in [E]J £ei a < (3 < jj, < 2a, 



\v\ a ir(dv) + / \v\^it{dv) < oo, 
v|<l J\v\>l 

and a\ tfi,g, f G C^(R d ),G^ G C* /3vl (R a! ). Iften /or each e G (0, 1) there is 
j3 > 2a and a unique u = u e G C^{H) solving (EI 



We separate in the operator L £ its "bounded jump" part L £ v(x) = L £ z v(x) 
with 



L%v{x) = -e a (-A) a/2 u + {a(z),V x v{x)) 
+\ Y,{b\z),V(z))dlv{x) 



2 



+ / [v(x + G (z) v) - v(x) - Xa(v)(Vv(x),G(z)v)] dn, 

J\v\<l 

z,x G R d ,v G C$°(R d ), so that 

L|v(a;) = L%v(x) + / [«(s + G(z)y) - v(x)]dir, x, z G K d . 

J\v\>l 

Remark 4. If the assumptions of Lemma [5| /ioM and ii £ G C^(H) solves 
\3. 3\) with (3 > 2a, then u £ satisfies the following equation as well: 

(3.4) (d t + L £ )u(t,x) = F(u,t,x), 

u[T,x) = g £ (x), 

where F(u,t,x) = F z (u,t,x)\ z=x with 

F z (u,t,x) = f(t,x) — / [u(t, x + G(z)v) — u(t, x)] dir. 

J\v\>l 
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Using a probabilistic form of a maximum principle we will derive uniform 
(independent of e) (7^-norm estimates of u £ and passing to the limit as 
£->0we will obtain u € C" (H) solving (|3.ip . First we prove some auxiliary 
statements. 

Let 

Z t = / / [(1 -x a {v))vp{dt,dv) + x a iv)vq{dt,dv)\ 

JO J\v\<l 

(3.5) = f f vq{dt,dv) + t [ (1 - Xa (v))vdn. 

JO J\v\<l J\v\<l 

For (s,x) G H,h G R ',^ G R rf , the following stochastic processes in [s,T] 
are used to derive the uniform estimates: 

dU t = edZf + a{Ut)dt + b{U t )dW t + G{Ut-)dZ t , 

dH t = [a(U t + H t ) - a(U t )]dt + [b(U t + H t ) - b{U t ))dW t 

(3.6) + [G(U t - + flt_) - G(C/ t _)] dZ t , 

^ = oW(C^ + fr t ;^)dt + 6W(^ + fr t ;y t )dWi 
+ / G< 1 >(U t _ + iZ t _;y t _)dZ tj 

7|w[<l 

d^ t = a (1 \u t -V t )dt + b^{UuVt)dW t + G^{U t ^V^)dZ u 
U s = x,H s = h,V s = Z,V s = Z, 

where Z a is R d - valued spherically symmetric a-stable process correspond- 
ing to (—A)"/ 2 and independent of Z. Recall for a function v on R rf we 
denote v^'(x;£) = (Vv(x),£),x, £ G R^ and, for example, componentwise, 

n m 

eft? = (Vai(U t ), V t )dt + £(V&*(l7 t )> V^cW? + ^(VG^(£/ t _), V t .)dZ l t , 

t=l i=l 

j = 1, . . . ,d. 

Lemma 3. (o^ If a 1 , &*• , G* J G C' 1 (R d ), i/ien /or each I > 2 i/iere is a 
constant C such that 

E[ sup |fr t |'] < C\h\ l . 

s<t<T 

(b) Ifa\b),G ij G C 1+K (R d ) with « G (0,1], then for each I > 2 £/iere is 
a constant C such that 

e[ sup |y s |'+ sup |y s | z ] < c|<e[ z , 

s<t<T s<t<T 

E[ sup l^-Ftl'] < CKI'M'". 

s<i<T 
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Proof, (a) Since ()3.5p holds, we have by Holder inequality and martingale 
moment estimates (see [IT], [T7] ) 

/rt \ V 2 r-t 

E sup \H r \ l < C[\h\ l + E[[ \H r \ 2 dr) ] + E[ / |# r | Z cir] 

< C[\h\ l +~E[[ sup |fl r /| i dr],a <t<T. 

Js s<r'<r 

and inequality follows by Gronwall lemma. 

(b) Similarly, for each I > 2, there is a constant C so that 

E[ sup \V s \ l + sup \V a \ l ] < C\£\ l . 

s<t<T s<t<T 

Then 

(ft \ V 2 ft 

/ |if r | 2K |y r | 2 (ir ) ]+E[/ \H r \ Kl \V r \ l dr] 
Js J Js 

(ft \ V 2 /-t 

y |F r - V r \ 2 dr\ ] + E[y |F r - J£| l dr] 

< C[E / \H r \ Kl \V r \ l dr + E f \V r — V r \ l dr],s <t<T. 

J s J s 

By Gronwall lemma, 

fT 

E sup \V r -V r \ l < CE / |iT r |^|K| z dr 

s<r<T Js 

< C / T [E(|^ r | M )] 1/2 [E(|K.| 2i ) 1/2 ]dr 

J s 

< C\C\ l \h\ Kl . 

□ 

3.1. Proof of Theorem |4j ^. Existence. By Lemma[2j for each e E (0, 1) 
there is a unique solution u e e C@(H) to flES) for some /3 > 2a. By Remark 
HI (|3.4p holds as well. Let (s,x) € if and C/t solves (|3.6p . By ltd formula, 

Eg £ (U T ) - u £ (s, x) = E f F(u e ,r, U r )dr 

J s 

and 

\u £ (s, -)|o < bio + y (|/(r, -)|o + C\u e (r, -)\ dr. 

By Gronwall lemma, there is a constant not depending on u £ and £ such 
that 

sup MV)|o < C[\g\ + / |/(r, -)|o(ir]. 

0<t<T Jo 
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As suggested in [7] , we estimate multilinear forms associated to the deriva- 
tives of u. Let k = [P]-, (t,x) G H,^,...,^ k G R d and 

il,..,i fe =l OX ^--- x H 

4°\t,x) = n £ (t,x). 
For z G R d , (i, x) e H^ 1 eR d ,...,£ k e K d , let 

P 2 4 fc )(t, x; e, ■ ■ ■ , £*) = -^(-A^u^ (t, x, £\ . . . , £ k ) 

+ / {«<*) (x + G(z)v; e + G« £>, . . . , £ k + G^ 1 ) £*» - 4 fc )(x; £\ . . . , £ fc ) 
^|«|<l 

-X»[(V,4 fc) (z; £1, • • • , G(z)v) - E(V 5 m( fc )(x; fr, . . . , £ fc ), G« (*; £»]}d7r 

z=i 

k 

+(a( Z ),v x 4 k \x; e,..., e)) + Y,( v e u(k) to* 1 * ■ ■ ■ ' £ fc )> ^ ^)) 

1=1 

+^{(6 i (z),^(z))44 fc )(x;C 1 ,...,^) 



+ E[(^( 1 )(z;^),^(z))^ui fc )(x; ,£\ ...,£ k ) 



1=1 



+ (b\z),^(z;e))d, lU ^(x;,e,---,e)}}- 

Differentiating both sides of (|3.4p and multiplying by ^ . . . we see 
that ui k \t, x; . . . , satisfies the equation 

(3.7) (t , x, e 1 , . . . , e fc ) + v s uf (t, x ,e,...,z k ) 

= A(u e ,t,x,^,...,^ k ), 

where 

A(u £ , t,x,e,---,H h ) = B(u £ , t,x,e,...,H k ) + F^(u e ,t, x; £\ . . . 
and B(u £ ,t, x, . . . , is a finite sum of the terms of the form 

[V^« (t, x + G(or)u; f 1 , • • • , C 1 ) - V x u® (i, x + G(x>; f S . . . , f ' )] x 

xG( w '(i;f«+ 1 ,...,^)!; 
-l 

d 2 u(p (t, x + sG(x)v; C 1 , • • • , f ' )G{x)vdsG {k - l) (x; f i+1 , . . . , f k )v 



with I < k — 2 and 

(t, x + G(x)i;; f 1 , . . . , f 1 )G Pl) (x; f V . . , ^ ) . . . G ( ' m) (x; ^ , . . . , ) 
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with m > 2, 1 < k, I + 1% + . . . + l m = k and (t? 1 , . . . , <f ; , . . . , (?%n ) being a 
permutation of . . . , In any case, there is a constant C independent of 
e and u £ so that for all (t,x) G [0, T] x R d ,f e R d , 

(3.8) lA^t.x.c 1 ,...,^)! < ccM^oifc + l/^OWIC 1 !...!^!, 

iA(u £ ,t, •,^ i ,...,^ fc )i/3-fc < c([« e (t,.)i^ + i/(t,oij8)i^[...ie fc i, 

and 

(3.9) \A(u B ,t, x,l\...,t)- A(u e ,t, x,e,..., H k )\ 

k 

< c(\f(t, ou + m*, oio E i^i • • • i^ii? - z l wz l+1 \ ■ ■ • l?D- 

On the other hand, for any (s, x) 6 H with the processes defined in (|3.6|) . 
it follows by ltd formula, 

B[4 k \r, u T , v^..., v$) - 4 k \ s , x, e, ■ ■ ■ , e fc )] 
= E[gi k \u T , v^..., v k ) - 4 k Hs, x, e, ■ ■ ■ , £*)] 

= vj\d t uf\t, U t , V t \..., V k ) + Vfj t 4 k \t, U t , V t \..., V k )]dt 
= B £ A( Ue ,t,U t ,V t \...,V k )]dt 

and 

E[4 fc )(T, U T + H T ,Vl..., V k ) - u«(T, U T ,v£,..., V k )} 

-[4 k \ s ,x+ h,e, ■ ■ ■ ,e) - 4 h) (s,x,e, . . . 

= E[g^(U T +H T ,V T ,..., V$) - gto(U T , V^,..., V$)] 

-i4 k \ s , x +h,e,---,e)- 4 k \ a ,x,e, . . 

= E / {[d t 4 k \t, U t + H t ,V t \..., V k ) + V e Ut+Ht 4 k \t, U t + H t , V t \ V k )] 

J s 

-[d t 4 k) (t, U t , V t \ ...,V k ) + Vfj t uf\t, U t , V t \..., V k )]}dt 
= E f [A{u e ,t, U t + H t ,V t \..., V k ) - A{u e , t, U t , V t \..., V k )]dt 

J s 

Since by fl£HJ| 

\A(u E ,t, U t , V t \..., V k )\ < C(\u E (t, -)\ k + \f(t, -)\ k )\V?\. . . \V k \, 
it follows by Lemma [3] and Holder inequality, 

(3.10) E\A(u £ , t, U t , V t \..., V k )\ < C{\u £ (t, -)\ k + \f(t, -)\ k )\e\ ■ ■ ■ \H k \. 
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\A(u £ ,t, U t + H t ,V t 1 ,..., V k ) - A(u £ ,t, U t , V t \..., V t k )\ 



< \A(u e ,t, U t + H t ,V t \..., V k ) - A(u E ,t, U t , Vf, V k )\ 

+\A(u e ,t, U t , V t l , V k ) - A(u £ ,t, U t , V t \ ...,V k )\ 
= Ai + A 2 , 

it follows by the estimates (|3.8j) , (|3.9() and Lemma [3] that 

EAi < C7(E|i7 i | 2 ^)) 1 / 2 (|/(t,.)| /3 + Kt,-)| /3 ) 
< C\h\P- k {\f{t,-)\ p + \u{t,-)\ p ). 

and for \h\ < 1 

k 

va 2 < c{\f{t,.)\ k + \u(t,-)\ k )Y,nv t l \--Avt 1 \\v}-vl\\vl +1 \... 

1=1 

< c(\f(t, -)\ k + \u(t, -)\ k ) Y,(Pl\v t l - ^l 2 ]) 172 !? 1 ! • • • I^ll^ 1 

1 

< C(\f(t,-)\ k + \u(t,-)\ k )\e\---\Z k \\hf- k . 
Similarly, we estimate 

n g f\u T ,v±,...,v k )\<c\g\ k \e\...\e\ 

and for \h\ < 1 

E\g( £ k \U T + H T ,V},..., V k ) - g( k \U T , V±,...,V k )\ 

< c\ g \ p \h\^ k \e\...\e\- 

So, 

\4 k \ s ,x-e,...,e)\o 

< C\S 1 \...\?\[\g\ k + J (\u e (t,-)\ k + \f(t,-)\ k )dt],0<s<T, 
and by Gronwall lemma, 

sup |4 fe )( s ,x;^ 1 ,...,^)|o<C|C 1 |...|^|[|gU+ [ T \f(t,-)\ k )dt}. 

0<s<T Jo 

Also, for \h\ < l,x € R d ,0 < s < T, 

\4 k \ s ,x + h,e, . . . ,e) - 4 k \ s ,x,e,. . . 

< c\hf- k \e\...\Z k \l\gb+ [ (\f(t,-)\p + Ht,-)\p)dt], 
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and by Gronwall lemma, 



sup \u (k \s,-^\...,C k )\fS-k 

0<s<T 

< c\e\...\H k \l\g\e+ [ \f(t,-)\pdt] 

Jo 



Therefore for each j3 G (a, 2a], 

(3.11) sup \u e \p<C[\g\ p + [ T \f(t,-)\ p dt]. 
ee(o,i) Jo 

Since for each (s, x) G H, 

(3.12) u £ (s,x) =g e (x)+J [L £ u e (t,x)- f(t,x)]dt, 

and there is a constant C > so that for all (i, x) G H,h G R d , 

(3.13) |^« e (i,x + /i) -d t u s (t,x)\ 

< \L £ x+h u(t,x + h) - L £ u £ u(t,x)\ + |/(t,s + fc) -f(t,x)\ 

< C\hf-<*(\u E \~ + \f\p) 

for some $ G (a, a + a A 1). It follows from (|3.1ip and (|3.13p that there is 
a sequence e n — > and w € C^(H) such that such K £n — » u uniformly on 
compact sets of i?. By (j3.11j) . L £ u £ (t,x) — » Lu(t,x) pointwise and passing 
to the limit in (I3.12p . we see that u G (^(if) is a solution to fj3. If) . 

2. Uniqueness. Let u 1 , u 2 G C^{H) be two solutions to (|3.ip . Then 
v = u 1 — u 2 satisfies (|3.ip with g = 0, / = 0. Let Xp x be the solution to 
(jl.ip starting from cc G R d at time moment s. Then by Ito formula, 

—v(s,x) = ~Ev(T, Xj, x ) — v (s, x) 

= E ( T [d t v(r, Xp x ) + Lv(r, Xp x )} dr = 

J s 

and uniqueness follows. 

4. One-Step Estimate and Proof of Main Results 



First, we modify the mollified function estimates for the Lipshitz spaces. 
Let w G Cg°(R d ), be a nonnegative smooth function with support in {\x\ < 
1} such that w(x) = w(\x\), x G R rf , and f w(x)dx = 1. Due to the symme- 
try, 

(4.1) / x i w(x)dx = 0, % = 1, . . . , d. 
Ju, d 

For x G R d and e G (0, 1), define w £ (x) = e~ d w (|) and the convolution 

(4.2) f £ (x) = J f{y)w £ {x - y)dy = j ' f(x- y)w £ (y)dy, x G R d . 
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Lemma 4. Let a < (3 < 2a, f G C^~ a (R d ). Then 

(4.3) \f(x) - f(x)\ < CeP- a \f\p_ a ,xe K d , 
and there is a constant C such that 

(4.4) \Lf\ < Ce^~ 2a \f\^ a 
Proof. Indeed, if /? — a < 1, then 

\f £ {x)-f{x)\ < J\f(x-y)-f(x)\vf(y)dy, 

< c\f\^- a . 

If p-ae (1,2], then 

\f(x)-f(x)\ = | / f{x + y) - f{x) - {Vf{x),y)w £ {y)dy\ 



< J J \{Vf{x + sy)-Vf{x),y)\dsw e {y)dydy 

< Ce?- a \f\^ a . 

According to Lemma 17 (iii) and Corollary 18 in |12| . for each /3 so that 

/3 — a < a 

\Lf\ < Ce^- a \f\^ a = Ce^ 2a \f\^ a . 

The inequality (14. 4p still holds for /3 — a = a or (3 = 2a by a straightforward 
estimate. □ 

We modify one-step estimate in |12j for Lipschitz spaces as well. 

Lemma 5. Let a < j3 < ji < 2a, 



\v\ a dir + / \v\^dir < oo, 

|i>|<1 J\v\>l 

and a\b^ G C^(R d ),G ij G C^ vl (R d ). Then there exists a constant C such 
that for all f G C p - a (R d ), 

\B[f(Y s ) - f(Y T J\T Tis ]\ < C\f\p. a 6i-\Sf8 G [0,T], 

where i s = i if ' Ti < s < t,, + i . 

Proof. Applying Ito's formula, for s G [0, T], 

E[f(Y s )-f(Y T J\T T J=E[ f {Ly Ti J%Y r ))dr\T Tis ]. 

Jr is 

Hence, for e G (0, 1), by (TO]) and (Oil . 

\W(Y s )-f(Y Tis )\T Tia }\ < |E[(/ - t)(Ys) ~ (/ - f £ )(Y T J\T Tia }\ 

+\B[f(Y s )-r(Y T J\F Tis ]\ 
< CF(e,5)\f\^ a , 
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with a constant C independent of e, / and F(e, 5) = e^~ a + e /3 ~ 2a S. Mini- 
mizing F(e,S) in e € (0, 1), we obtain 

\E[f(Y s ) - f(Y T J\T T J\ < C^liV 

□ 



4.1. Proof of Theorem [H Let u € &{H) be the unique solution to (|3.1 
with / = 0. By ito's formula, 

EK0,X )] = B[u(T,X T )]--E[ [ T (d t u(s,X s ) + L Xs u(s,X s ))ds] 

Jo 

= E[<7(X T )] 

and 

(4.5) EK0,Xo)] = E[u(0 } y o )]- 
By Lemma [U 

(4.6) |L z n(s,-)|/3- a < C\ g y,\d t u(s,-)\p- a < C\g\p,s£ [0, T], z £ K d . 
Then, by Ito's formula and (|4.6p . it follows that 

B[g(Y T )} - E[g(X T )] = E[u(T,Y T )] - E[u(0,y )] 

= E 



{[d t u(s,Y s ) - d t u(s,Y T J] 
+ [Ly Tis u(s, Y s ) - L Ytis u(s, Y Ti3 )] }ds 



Hence, by (|4.6p and Lemma [5l there exists a constant C independent of 
g such that 

\Eg{Y T )-Eg(X T )\<C5^- l \g\ p . 
The statement of Theorem [1] follows. 

4.1.1. Proof of Corollary^ According to [2J, there is a rapidly decreasing 
smooth function w € <S(R d ), the Schwartz space, such that f w(x)dx = 1 
and all moments are zero: 

where x 7 = xj 1 . . . x^ d , x = (x\, . . . ,Xd) E R d . Let e E (0, l),w £ (x) = 
e~ d w(x/e), x G R d , 

& (^) = j g(x- y)w £ (y)dy, x e n d . 

We will show that for G (0, 4] , v < 0, 

(4.7) sup\g £ (x)-g(x)\ < C\g\ v e\ 

X 

\g e \p < Ce v ~%\ v 
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(A standard mollifier could be taken if v < 2, see Lemma S]). Since for 
x G R d , 

9e{x) - g{x) = I \g{x - y) - g(x) - ^ — 9 ^ y J ]w £ (y)dy, 
it follows that 

sup \g £ (x) - g{x)\ < C\g\ v e v . 

X 

If /3 is an integer, 7 G N d , [ — y | = /3 and 7 = /x + // with |/i| = [1/], // 7^ 0, 
then 

D^ e (x) = j g(y)(D^w) £ {x-y)dy = e^ j D»g{y){D»' W ) £ (x-y)dy 

/" [^(y) _ ^^^(a;)] (D^'w) E (x - y)dydy 



and 

|£> 7 g e (x)| <^|i,x6R d . 

If /3 is not an integer, the second inequality in (|4.7p follows by interpolation. 
According to Theorem Q] and (|4.7p , 

\Eg(Y T ) - Bg(X T )\ < 2 sup \g £ {x) - g{x)\ + |E 5e (y T ) - Eg £ (X T )\ 

X 

< C\g\ u F(e,5), 

where F(e,5) = e u + e u ~ /3 5^~ l . Minimizing F in e G (0, 1), the statement 
of Corollary [2] follows. 

4.2. Approximate simple Euler scheme. Consider the approximation 
of Xt defined by the increments of Zt = Z" + R^, < t < T, in Example [TJ 
Obviously, Z t depends on a, /3 and a. Its generator is 

Lv(x) = / / a; + i;) — f (s, x) — x Q (f ) (Vf (s, x), u)]7r((if ) 

•/>!>£ 

+ J R a ' /3 u(x), 

where 

f ^(V^),^ if a < /3 G (1, 2], a G (0, 1], 

if a </? G (2, 4], a G (1,2], 

otherwise. 

Lemma 6. Let a < (3 < 2a and h G C^(R d ). Then there is a constant C 
such that for every ¥ za -stopping times < r < r' < T we have 

\E[h(Z T > - Z T ) - h(Z T , - Z T )\F T ]\ < C0(o-)|%E[r' - t\F t ], 



V[X 
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with 

0(<r)= / \vf A3 dn 



iiKu 



(Ttere F^ CT is the natural filtration of a -algebras generated by Z a ). 

Proof. Let Z u = Z — Z° ' . We show first that there is a constant C such that 
for any s < t,g G C p (R d ), 

(4.8) |E 5 (Zf - Z%) - Eg(R? - < C<f>(a)\g\ p \t - s\. 
By Ito formula 

(4.9) v(r,x) = E\g(Z? - Z° + x),0 < r < t, 
is the solution of the backward Kolmogorov equation 

(4.10) dtv(r,x)+ / [v(r, x + v ) — v(r, x) — X a (. v ) (Vw(r, a;), v)]ir(dv) 

J \v\<<7 

= 0,v(t,x) = g(x),0 < s < t. 

Obviously, v G C^([0,t] x R d ) and (see (|Q]> ) \v\p < \g\p. By Ito formula 
and dUD]) , 

(4.mRt-Rs)-Vg(Z?-Z°) 
= Bv(t, R1 - R a s ) - v(s, 0) = E J [R a ^v(r, R° - R°) - Lv(r, R u r - R°)]dr, 
where 

Lv(r,x) = / [v(r, x + v) — v(r, x) — x a ( v ) (Vi>(r, x), v )]ir(dv), (r, x) G H. 

J \v\ <cr 

If a < P G (1,2], a G (0, 1], then for all (r,x) G H, 

\R a '^v(r,x)— / [v(r, x + v) — v(r, x)]ir(dv)\ 

J \v\ <cr 

< / / |V«(r, x + st> ) — Vv(r, x)\ \v\dirds 

JO J\v\<a 

< C\v\ p [ \vfdir < C\h\ p [ \vfdir. 



V\<(7 J\v\<(T 



If a < P G (2, 4], a G (1,2], then for all (r,x) G H, 

\R a '^v{r, x) — I [v(r, x + v) — v(r, x) — (Vv(r, x), v)]dTr\ 

J\v\<a 

n\D 2 v (r, x + sv ) — D 2 v(r, x)\ \v\ 2 dnds 
>\<a 



< 



< C\v\ p / \vf A3 d7r < C\h\ fi / \vf A3 dn. 

J\v\<a J\v\<a 
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The estimate of the difference R a '^v—Lv in the other cases is straightforward 
and @3D follows by grp. 

Since Z a ,Z a and R a are independent and r, r' are F^ CT stopping times, 
we have by $M that 



\E[h(Z2 -Z° + Z%, - Z°) - h(Z°, -Z° + R%, - R°)\T T ]\ 
< C<t>(a)\h\{,E[T J -T\T T \. 
The statement follows. □ 

For the proof or Theorem [2] we will need the following estimate. 
Lemma 7. Let 

Vt = at + bW t + GZ t , 

where a G R d , b is a d x d-matrix and G is a m x m-matrix. We assume 
b = if a G (0, 2) and a = if a E (0, 1) and 

|o| + |6| + |G| <if. 

Leta< /3 <fi<2a andhe C^~ a (R d ). 

T/ien i/iere is a constant C = C(a, f3, K) such that 

\Bh(V t ) - h(0)\ <c£- l \h\ p _ a . 

Proof. For / 6 (7^(R d ), applying Ito formula, 

E/(Vt) - /(0) = E f Kf(V r )dr, 
Jo 

where for x G R d , 

/C/(x) = (a,V/(x)) + ^6*6^./(x) 



+ J [f{x + v)- f{x) - xM(Vf{x),v)Wv). 

For h G 6P- a (B. d ) we take u; G Cg°(R ), be a nonnegative smooth function 
with support in {|xj < 1} such that w{x) = u>(|x|), x G R d , and f wix)dx = 
1. For x G R rf and e G (0, 1), define u> e (x) = (|) and the convolution 

h £ ix) = J fiy)w £ ix-y)dy,x G R d . 

Then by Lemma H] 

\Eh(V t ) - h(0)\ < 2e p ~ a \h\p_ a + |E / O 

< C|%_ a (^- a + e p - %a t) 

for each e G (0, 1). The statement follows by minimizing the inequality in 
e. □ 



t 

£ {Vr)dr\ 
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4.2.1. Proof of Theorem® Let u G C@(H) be the unique solution to the 
backward Kolmogorov equation 



(4.12) 



Let for Tj < t < Ti+i 



(d t + L)u(t,x) = 0, 

u(T,x) = g{x). 



HI = aiY^it-T^ + biY^Wt-Wr,) 
+G(Y n )(Z t -Z Ti ) 

and denote AY Ti = Y Ti+1 — Y Ti . We approximate 

u(T,Y T ) - u(0,Y ) 
= /~2u(T i+1 ,Y Ti+x ) -u(ri,Y Ti ) 

i 

= J2Hr i+ i,Y Ti + AY Ti ) - u(r i+1 ,Y Ti + H^ +1 )] 

i 

+ J2Hr i+1 ,Y Ti + lP Ti+1 ) - u(n,Y Ti )} 

i 

= D 1 + Y J D 2i - 

i 

According to (12.4j) (Lemma [6} , 

E|Di| <C<t>{a)\u\p <C<K<r)\g\p. 
Now, we estimate the second term. By Ito formula for each i, 
E[D 2i \T Ti ] = B[u(T l+u Y Tt + Hl t+i )-u(r i+u Y Tt )\T Tt } 

= E{ r + \d t u{rX i + Hi)+Ly r u{r,% i + H i r )]dr\F Ti } 

= E l(d t u(r, Y Ti + H l r ) - d t u(r, Y Ti )) 

■JTi 

+ {Ly T u(r,Y Ti + K) - L ?t u(r,Y Tl Wr 
and by Theorem [H and Lemmas [I] and 



< J^|ED2i| < C8«- x \Lu\p. 



and the statement of Theorem [2] follows. 



2-1 
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4.3. Approximate jump-adapted scheme. Consider the approximation 
of Xt defined by the increments of Zt = Z° + Rf,0 < t < T, in Example 
CD For a £ (0,1), 5 > 0, consider the following Z CT -jump adapted time 
discretization: tq = 0, 

T i+1 = inf (t > n : AZ? + 0) A (t< + 5) A T. 

In this case the time discretization {Tj,i = 0, . . . ,ny} of the interval [0,T] 
is random, are stopping times. We approximate Xt by 

y t = X + I a(Y T Jds + [ b(Y T JdW s + f G(Y T JdZ s ,t € [0,T]. 
io io io 

In this case, 

Ti+i - Ti = r/ i+1 A 5 A (T - Ti) 

with 

=mi(t>0:p((Ti,Ti + t],{\v\ > a}) > 1) 
and r] i+1 is J 7 -^ -conditionally exponential with parameter \ a = ir ({\v\ > a}). 

Lemma 8. Let 5^ = 5 A (T — Tj), i > 0, and Ao- = 7r > a}). 
(%) There is constant c > suc/i i/iai /or any i > 

c $ A A; 1 ) < E[r m - n \T Ti ) <5' t A \-\ 

(ii) There is a constant C such that for any i > 0, 

Ef^+x-Ti) 2 !^] < CB[5[ 2 A\- 2 \T Ti ] 

< C , ($AA^ 1 )E[r i+ i-T i |.F T ,]. 

Proof. Since Tj+i — Tj = A 5 A (T — Tj) and 

r) i+1 = inf(i > :p((Ti,Ti + t],{\v\ > a}) > 1) 

is T Ti -conditionally exponential with parameter \ a , we find 

r s 'i 

E[T i+ i-Ti\T Ti ] = E[ Vi+1 A8i\T Ti ] = \ a / te'^dt + ^e-^ 

Jo 

1 — e -A CT <5^ 

a; • 

If ^ > A" 1 , then 6% > 1 and 



1 - e _ACT,5 « 1 - e" 1 1 



> — > -A 



-i 



A(j Ao- 3 

If 5- < A" 1 , then ^A CT < 1 and 
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Therefore (i) follows. Similarly, 

fS'i 

E[(r i+ i - T t f\T Ti ] = X a E[ Pe-^dt + Sfe-^lTrJdt 

Jo 

-V 

and (ii) follows using (i). □ 

An immediate consequence of Lemma [8] is the following statement. 
Corollary 5. (i) There are constants c, C > such that 

cEj> i+ i-Ti) < A A" 1 ) A {T-Ti)} 

i i 

< CE^(r !+1 - n) = CT. 

i 

(ii) There is C > such that 

^E[(r m -r,) 2 ]<Cr(<JAA^). 

i 

Proof. We derive (i) by summing inequalities in Lemma EJi). According to 
Lemma E^ii) and (i), 

^E[(r m -r 4 ) 2 ] < C^E[(T-r,) 2 A£ 2 AA; 2 ] 

i % 

< C(TA5A\- l )Y,n(T-Ti)A5A\- 1 } 

% 

< CT{8 A A" 1 ). 

The statement follows. □ 

For the proof of Theorem [3] we will need the following estimate as well. 
Lemma 9. Let 

V t = at + bW t + GZ t , 

where a G R d , b is a d x d-matrix and G is a m x m-matrix. We assume 
b = if a £ (0,2) and a = i/ a E (0, 1) and 

M + |b| + \G\ < K. 

Leta< (3 < /x < 2a and h £ C^ _a (R d ). 

TTien i/iere is a constant C = C(a, f3, K) such that for any i > 

|E[ ^(^-WJI^rJI < C\h\^ a \f -1 [5 A X- 1 )-- 1 E[(r i+1 -r 4 )|J- T J, 

J Ti 

where \ a = ir ({\v\ > a}) , 

X a = 1 + l ae (i i2 )| / vdn\. 

Jl>\v\>a 
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Proof. For / € C^(R d ), i > 0, applying Ito formula, 
E 



[ f(V r )-f(V Ti )\T Ti ]dr = E [f Kf(V r )dr+M s -M Tt ]ds\F Ti ]dr, 

where for x G H d , 

Kf{x) = (a,V/(x)) + ^6*6^/(x) 

+ J [f(x + v)- f(x) - Xa (v)(Vf(x),v)]ir(dv) 



and 



Mt= 11 I [f( - Vr ~ + Gv) ~ f( - Vr -^ q ( dr ' dv ^ f e [°> ^ 



Note that 



/ (M s — M Ti )d{s — Tj) = (M ri+1 - M r J(r j+i — Tj) — / -Ti)dM a . 

Since Z a &nd Z a = Z — Z a are independent and Tj are -stopping times, 
it follows by definition of Tj that 

E[(M ri+1 - M Ti )(r i+ i -n)- f (s - Ti)dM s \F Ti ) 

J n 

= E[-(T i+1 -T i )(U^ +1 -U^)+ r + \s-Ti)dU^ Ti ] 
= -E[ r + \u°-U° Ti )ds\F Ti l 

J Ti 



where 

<-t 



U? = [ [ [f(V r - + Gv)-f(V r -)]dTrdr 

J J\v\>cr 

= f I [f(V r - + Gv) - f(V r -)]dTrdr + f I x a (v)(Vf(V r ),v)d7rdr 

JO J\v\>l JO Jl>\v\>a 

+ f l [f(V r - + Gv) - f(V r -) ~ X a (v)(Vf(V r ),v)]dTTdr 

JO Jl>\v\><7 



Hence 



(4-13) |E[ f T +1 f(V r ) - f(V Ti )\T Ti \dr\ 

J Ti 

< C(l + l ae(1)2) | / vdn\)\f\ P mT i+1 -Ti) 2 \T Ti ]. 

Jl>\v\>e 

For h £ C /3 ~ a (R d ) we take w € Cg°(R d ), be a nonnegative smooth func- 
tion with support in {\x\ < 1} such that w(x) = w(\x\), x € R d , and 
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f w(x)dx = 1. For x G H d and e G (0, 1), define w e {x) = e d w (§) and the 
convolution 

h £ (x) = J f(y)w e (x-y)dy,xeK d . 
Then by Lemma 0] and (|4.13p , 



E[ 



h(V r ) - h(V T )\F Ti ]dr 



T i+1 



< 2£^- a \h\^ a B[(r i+1 -T i )\F T ] + \E[ I (h £ (V r ) - h £ (V T ))dr\T T }\ 

J n 

< 2^- Q |%_ a E[r i+1 - Ti\T Ti ] 
+Ce^ 2a [l + l ae{1 , 2) \ [ vd-Km^niri+i-rif^). 

Jl>\v\>e 

Minimizing the inequality in e we find by Lemma On) that 



E[ 



Ti+l 



h{V r ) - h{V T )\F Ti ]dr 



< c\h\^ a x; 'e^+i - Ti\T T f-^[(ji+i - Tif\T Ti ]«- 1 

< C\h\ p _ a K~ l (8 A A^^Efo+i - Ti\F Ti ]. 



□ 



4.3.1. Proof of Theorem^ Let u £ &(H) be the unique solution to the 
backward Kolmogorov equation (see Theorem U|) 

(4.14) (d t + L)u(t,x) = 0, 

u(T,x) = g(x). 

Let for Tj < t < Tj + i 

+G(Y Ti ){Z t -Z Ti ) 
and denote AY Tt = Y Ti+1 — Y Ti . We approximate 

u(T,Y T )-u(0,X ) 
= J2u(T i+1 ,Y Tx+1 ) - u(n,Yr.) 

i 

= J>(r i+ i,r T4 + AY Ti ) - u(r i+1 ,Y Ti + K l+1 )\ 

i 

+ Y Ti + K i+1 ) - u(r u Y Ti )} 

i 

= D 1 +^ j D 2i . 
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According to Lemma [U 

E|Di| <C<l>{a)\u\p<C<t>{a)\g\p. 
Now, we estimate the second term. By Ito formula for each i, 
E[D 2i \T Ti ] = B[ U (T t+1 ,Y Tt +K t+i )-u(T t+1 ,Y Tl )\^} 

= E{ [d t u(r,Y Tl + K) + L f u(r,Y Ti +Hi)}dr\T Ti } 

= E l(d t u(r, Y T . + HI) - d t u(r, Y Ti )) 

J Ti 

+{L % u(r,Y Ti + K) - Ly u(r,Y T J)]dr 
and by Theorem [J] and Lemmas [TJ [9] and Corollary [5j 



^EL» 2i 



< ^2\ED2i\<CXS '(^AA; 1 )" 1 {\d t u\ p _ a + \Lu\p_ a ) 



< 



CAf ^SAX- 1 )' X \u\p<C'4 '{SAX- 1 )- l \g\ p 



and the statement of Theorem [3] follows. 

5. Conclusion 

The paper studies a simple weak Euler approximation of solutions to 
possibly completely degenerate stochastic differential equations driven by 
Levy processes. The dependence of the rate of convergence on the regularity 
of coefficients and driving processes is investigated under the assumption of 
/3-Lipshitz continuity of the coefficients. It is assumed that the SDE is 
driven by Levy processes of order a € (0, 2] and that the tail of the Levy 
measure of the driving process has a /i-order finite moment (/i G (a, 2a]). 
The resulting rate depends on f},a and \i. Following [3], the robustness of 
the results to the approximation of the law of the increments of the driving 
noise is studied as well. It is shown that time discretization and substitution 
errors add up. In addition, a jump-adapted approximate Euler scheme is 
considered as well. The derived error estimate shows that sometimes the 
inclusion of jump moments into time discretization {t{\ could improve the 
convergence rate. In order to estimate the rate of convergence, the existence 
of a unique solution to the corresponding backward degenerate Kolmogorov 
equation in Lipshitz space is first proved. 

On the other hand, there is a discrepancy in the model (II. ip between 
a = 2 and a € (0, 2). One would like to consider the equation 

X t = X + [ a{X s )ds+ [ b(X s )dW?+ [ G(X s -)dZ s ,t G [0, T], 
Jo Jo Jo 

with a possibly degenerate b and a spherically symmetric a-stable W a (in 
(fTTjl . 6 = for a € (0,2)). 
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Since (jl.ip could be degenerate, a solution corresponding to a given a G 
(0, 2] can be looked at as a solution corresponding to a € (a, 2] as well. 
Therefore the rate for a fixed a cannot be "universally optimal" : there is 
always a large subclass for which the rate claimed for a could be better 
and achieved under weaker assumptions. For example, if /3 = = 2a with 
a S (0,2) (the diffusion part is absent), the convergence order is re = 1 
(fi = 4 and G € C 4 is not needed). Even "strictly at a", the assumption 
about the tail moment fi G (a, 2a] is not optimal. It could be weakened for 
a subclass with the driving processes Z such that the compensator of the 
jump measure of Xt has a nice density with respect to a reference measure. 
For example, let us consider the following one dimensional model 

(5.1) X t = X + [ a{X s )ds+ [ b(X s )dW s + [ G(X s -)dZ s ,te [0,T], 
Jo Jo Jo 

where Z is a symmetric A-stable with A S (0, 1) and G > 0. Assume 
a,b,G x ,g € (7 4 (R). Although fj, < 1 in this case and the equation is pos- 
sibly degenerate, a plausible convergence rate is still re = 1 (or re = if 
g € C^(R), v € (0, 4]), because the integral part of the generator of (|5.1j) . 

/dy f dy 

[v(x + G(x)y) - v(x)] +A = G(x) x J [v(x + y) - v(x)\ i 

is differentiable without assuming much about the tail moments of the Levy 
measure. 
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